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^ ! Abstract 

£C) . Using tidal disruption of globular clusters by the galactic center, we 

' put limits on the total mass ever enclosed into orbits of observed globular 

t-H i clusters. Under the assumption that the rate of mass loss from the Galaxy 

■ is steady, we then deduce a bound on this rate. In particular this bound 

can be used to constrain the galactic gravitational wave luminosity. 

in ■ 

o : 

^ ; 1 Introduction 

i-S^ ■ Upper bounds on the mass loss of the Galaxy have been discussed in the past, 

especially in connection with the maximum luminosity of the Galaxy in gravi- 
tational waves [QI2]- Recently, in Ref. [3], the issue has been revisited, partly 

£j ; motivated by the 2001 data of the NAUTILUS/EXPLORER collaboration 0). 

As discussed in [3J, limits on the mass loss are given by its effect on galactic 
dynamics. If the Galaxy loses mass, it becomes less bound and expands. One 
can constrain this expansion from various observations; in particular: 

Jv>( , (i) The radial velocity of stars in the sun's neighborhood. In particular a com- 

5^ ■ ponent linear with the distance and independent of the direction (i.e. a K-term 

in the galactic dynamics literature). 

(ii) The relative velocity of this neighborhood with respect to the galactic cen- 
ter. 

(iii) The existence of old wide binaries near the galactic center. 

(iv) The existence of globular clusters near the galactic center. 

This last method, based on tidal disruption of globular clusters, was proposed 
by A. Poveda and C. Allen [5] and applied to the globular cluster Omega Cen- 
tauri. Actually some facts about its metalicity, mass and shape support the 
idea that oj cen is a rather unusual globular cluster and probably the remains of 
a tidally stripped dwarf galaxy |S[ [3 |Hj ■ In this case, it would not satisfy the 
assumption of an adiabatically evolving orbit required in the analysis of Ref. 
PJ. For this reason we repeat this analysis using other globular clusters. 
We first discuss in Sec[3 the evolution of elliptical orbit in the potential of a 
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varying mass. Then wc recall some facts about the tidal radius and the dynam- 
ical evolution of globular clusters. In Sec|3 combining this and the observations, 
we are able to derive an upper limit on the mass loss. 



2 Evolution of the globular clusters 

In this section we derive the tidal radius of globular clusters using Newtonian 
gravity and disregarding other effects. 

If there is a continuous mass loss in the Galaxy, stellar orbits change with 
time. However, for slow changes of the galactic mass, there exist some adiabatic 
invariants in the description of the orbit (see [5] Chap. 3). In the case of an 
ellipse around a mass M, the eccentricity e and the product of the mass with the 
major semi axis, Ma, are such invariants. This adiabatic invariance property 
can be then extended to more general cases ^U] ■ For an axially symmetric field 
and orbits on the galactic plane the adiabatic invariance of the eccentricity is 
still valid (with a more general definition for the eccentricity for non elliptical 
orbits). In particular, if we consider the osculatory ellipse at perigalacticon 1 , we 
have M = M p (where M p denotes the mass enclosed in a sphere of radius R p ), 
and the adiabatic invariance means e = constant and M p a — constant, from 
which follows also M P R P — constant, which is the formula we will use in the 
following. As mentioned above, this is correct if the motion lies in the galactic 
plane; in general this is not true for globular cluster orbits. However, if the 
orbit is of a box type there is on average no exchange of energy between the 
vertical motion and the motion in the plane the two motions can thus be 
separated and the invariance theorem can be applied to the motion in the plane. 
For this reason we consider in the following only globular clusters which have 
box orbits. A star cluster orbiting in a potential is subject to tidal forces. Due 
to these forces the cluster could lose its most distant stars. We follow King's 
definition of the tidal radius ^2] Let's consider a spherical globular cluster of 
mass M c describing an ellipse around the galactic center. The acceleration of 
the cluster at a distance R with respect to the galactic center is given by: 

d 2 R 2 dV 

where u> is the angular velocity and V the galactic potential. For a generic star 
in the cluster, the acceleration is: 



dt 2 V dR 



M C {R*-R) 
G \R-RJ3 ' (2) 



where R* is the position of the star measured from the galactic center. 

Using a Taylor expansion we can compute the relative acceleration to first order 

1 For non elliptical orbits the osculatory ellipse at perigalacticon is the ellipse approximating 
the orbit near this point. 
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in - R): 



d R* d 2 R ( 2 d 2 V GM i 

uj — 



i 

tfta dP ~ V~ di? 2 |A. - R\ 3 J {R * ~ R) ' (3) 

The tidal radius is defined as the distance from the cluster center for which the 
relative acceleration vanishes (it corresponds to the maximum distance that a 
star can reach without escaping from the cluster): 

rt -\^m ■ (4) 

Indicating with Mqc the galactic mass inside a radius R and assuming that it 
is spherically distributed we have: 

* V 2G Mgc (5) 

Using the geometry of the ellipse, ui 2 — GMccR ia {^ — e2 )i we obtain for the 
tidal radius at perigalacticon 2 (R = R p = a(l — e), Mqc = M p ): 



'•'=Ms^)J ■ (6) 

In the following we make use of a well known fact about globular clusters 
evolution: in the first phase of the evolution, until core collapse, the central part 
of a globular cluster undergoes a dynamical contraction (see for example |14j- 
|18p while the outer part expands. The globular clusters that we are considering 
in our analysis are all in this evolutionary phase. It remains difficult to determine 
the fraction of the mass that makes the contraction. Poveda and Allen |2] 
estimated it to be 75% making use of results found by Spitzer and al. |141 ITS] 
using numerical simulations. According to more recent estimations II 7| we 
consider the smaller fraction 1/2. Actually for our purpose the only thing we 
need to consider is a fraction of the cluster mass which undergoes a contraction. 
It is clear that if we consider the maximal fraction with this property we obtain 
the tightest bound on the mass loss; but in any case even with a smaller fraction 
we obtain a limit, even if is less strict. This explain why we obtain a different 
bound than Poveda and Allen for ween. 



3 Limit on the total mass loss 

The tidal radius of globular clusters gives informations on the mass of the galac- 
tic center enclosed in a sphere of radius R p . By mass we mean total gravita- 
tional mass and therefore the analysis makes no difference between "normal" 

2 The tidal radius at perigalacticon can be considered the effective tidal radius of the cluster: 
for a discussion of this issue see 1121 1131 
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and "dark" matter. Furthermore the limit on mass loss gives no contraint on 
the nature of the mass lost or on the loss mechanism. 

Following the idea of Poveda and Allen [2] , we can obtain a limit on the total 
mass loss, assuming at first that the mass loss is localized at the galactic center. 
As described above, the radius r x , containing a fraction x of the cluster mass, 
decreases with time for small enough x. We have the inequality: 

r\ = ar l x > ar™ , (7) 

where a is the quotient between the tidal radius r t and r xi the superscripts i 
and n stand for "initially" and "now" respectively. 

Using the formula for the tidal radius and the adiabatic invariance of M p R p and 
e we have: 

&l( — V /3 > ar » (8) 

From this equation we can derive an upper limit on M*. In order to do this, we 
should use a model for the mass distribution to determine a, but it is clear that 
it is in any case smaller than for the case of a uniform density distribution, for 
which a = l/x 1 / 3 . Thus we have: 



,1/3 ' X 



Since all the galaxy has expanded then M™, the mass enclosed in a sphere 
of radius R p , was originally enclosed in a sphere of radius R p and therefore the 
mass loss is given by: 

AM = AP p - M p l . (10) 



3.1 Results 

The most stringent limits correspond to clusters with small peri galacticon ra- 
dius. On the other hand clusters with R p < l[kpc] have chaotic orbits and do 
not satisfy the adiabatic invariance conditions. We restrict our choice to clusters 
known to have box orbits, see [T^l, but with R p small enough; in all the cases 
that we consider, R p is approximatively the size of the bulge, so we can take 
M p as the mass of the bulge. 

Orbit type determination are obtained taking a fixed galactic potential (this 
is also the case in It would be interesting to study how the result changes 

if one use a time- varying potential wich takes into acount the variation of the 
mass of the Galaxy. However, we do not expect substential changes in the final 
result because of the structure of equation iJSJ . In fact a miss-determination of 
R p and e of said 20% do not change the order of magnitude of our result. 

The parameters that we need for our selected clusters are: 
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Name 


NGC 


M c [Mq] (:20 ) 


R; [kpc] ( 21 ) 




ro.5 [pc] ( 22 ) 


uj cen 


5139 


2.40 x 10 e 


1.2 


0.67 


11.7 


M107 


6171 


6.31 x 10 4 


2.3 


0.21 


7.54 




6218 


1.17 x 10 b 


2.6 


0.34 


7.96 


M92 


6341 


2.19 x 10 b 


1.4 


0.76 


5.75 


M28 


6626 


2.29 x 10 b 


2.1 


0.19 


7.17 



With Mb = 1.6 x 10 10 [M Q ] the mass of the galactic bulge, we obtain 





ui cen 


M107 


NGC6218 


M92 


M28 


AM [Mq] < 


1.9 x 10 10 


1.7 x 10 10 


2.4 x 10 10 


2.0 x 10 10 


2.7 x 10 10 



If the mass is only lost from the center of the galaxy, then this is the total mass 
loss. For a steady loss during the whole lifetime of the galaxy (and taking the 
lifetime of the galaxy to be At = 1.2 x 10 10 [yr]) the maximum rate would be 





lu cen 


M107 


NGC6218 


M92 


M28 


^-[Me/yr] < 


1.6 


1.4 


2.0 


1.7 


2.3 



From the above analysis we see that only the variation of the mass inside the 
globular cluster orbits is taken into account, see eq.©. We have assumed that 
the mass loss takes place only there. For example the mass loss could be pro- 
portional to the luminous mass distribution. In this case, the total mass loss is 
given by 

AM tot = —AM (11) 
M B 

where M tot denotes the total luminous mass in the galaxy, AM the mass loss 
in the bulge previously calculated, and AM 4ot the total mass loss in this case. 
Taking M fot = 1.1 x 10 11 [M Q ] we obtain, assuming a steady loss during the 
whole lifetime of the galaxy, the limits on the rate of mass loss 





uj cen 


M107 


NGC6218 


M92 


M28 


M^{M G /yr] < 


11 


9.7 


14 


11 


16 



4 Conclusion 

Globular clusters allow to fix a limit on (a steady) galactic mass loss of order of 
l[M©/?/r] for the galactic center (central kpc) or, if we assume that the mass loss 
is correlated with luminous mass, of order of 10[Mo/yr] for the whole galaxy. 
These limits are independent of the kind of mass loss and then could be used 
to constrain any mechanism producing a net outflow of gravitational mass. 
In particular the first limit gives a bound on the maximal gravitational wave 
luminosity from the galactic center. In 0] has been reported the analysis of the 
data of the 2001 run of the EXPLORER and NAUTILUS gravitational wave 
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detector. More statistics is needed to draw any conclusion (and further runs are 
presently under way), but in that data there was a hint of a possible coincidences 
eccess when the bars are oriented favorably with respect to the galactic plane. 
If interpreted in term of gravitational wave emission from galactic sources, this 
would correspond to a mass loss of the Galaxy to gravitational waves of the 
order of (2 — 6) [M® / yr] , see |3J . Combined with other arguments (see [3] for a 
detailed discution) , our work shows that a model with sources of gravitational 
wave located in the galactic center, radiating steadliy for a time comparable 
with the age of the Galaxy, is only marginaly consistant with the observation. 
For this reason in [3] have been considered other alternatives, like a single or 
just a few close sources emiting repeatesdly gravitational wave burst. 
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